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Abstract 

The effect of confinement on surface plasmon polariton in a planar nanocavity 
was studied. Tlie generalized modes were obtained and studied in detail. 
It was demonstrated that these modes result from the strong coupling of 
plasmon-like and photon-like modes. 
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1. Introduction 

Surface Plasmon Polaritons are surface mode localized at the interface 
between a dielectric and a metal that can be confined in regions much smaller 
than their wavelength. Because of this strong enhancement many interesting 
phenomena are linked to plasmons: a huge Purcell effect [1], the strong 
coupling with excitons [2], quantum dots luminescence [3], spasers [4]. All 
this has resulted into the new field of Plasmonics [5]. It has even be suggested 
that plasmons could allow for a control of light at the quantum level [6] . 

In this work, we study the effect of confinement on surface plasmons in 
a planar cavity [7] with a wall coated with a lossy metal. The dispersion 
curves are studied in details and explicit relations are obtained locally. In 
particular, it is shown that the modes result from the strong coupling of 
photon-like modes with plasmon-like modes. It is shown that modes strongly 
concentrated at the air-metal interface can be obtained. Finally, we discuss 
the modes of the lossy waveguide, in particular their completeness. 
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Figure 1: The structure under study. 

2. The dispersion relation 

The structure under study is depicted in fig. (1). We consider H\\ modes 
only, that is electromagnetic fields for which the magnetic field is linearly 
polarized along the z direction. Time harmonic field with a time dependence 
of e"*^"'^* are considered (c is the speed of light in vacuum). We denote 
H{x,y) = u{x)e^'^y e~^^°'^* Cz- The total heigth of the waveguide is h and the 
thickness of the metal is d. We denote t = d/h. The field u{x) satisfies the 
following equation, in the Schwartz distributions meaning: 

(e-'d^u) + {kl - j'e-') u = (1) 

2 

where e{x) is equal to 1 for x G [d, h] and to Em = '^~zj(Zj+ir) ^'^^ ^ ^ l^^d], that 
is, the metal is described by a Drude model. For numerical computations, 
the parameters of silver will be used (the plasma wavelength is 137nm). As 
a simple model of confinement, we use the Neumann conditions dxU = on 
the boundaries of the cavity (x = 0,h). 

The plasma frequency Up allows to define a unit of length by considering 
the associated wavelength A = ^ and a unit of spatial frequency by con- 
sidering the corresponding wavenumber fcp = We denote with the upper 
case the normalized length quantities: H = kph the normalized height of the 
waveguide, D = kpd the normalized width of the metal layer, X = kpX the 
normalized variable, Kq = ko/kp, G = 'y/kp the normalized wavenumber and 
propagation constant respectively. We denote 
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and the corresponding dimensionless quantities: 

Pi = Kl -G\Pl = Kle^ - (2) 

A straightforward computation leads to the following expression for the 
field: 

Q < X < d : u{x) = B cos{(5rnX) 
d < X < h: u{x) = Ccos{l3g{X - H)) 

The transmission conditions at the interface X — D, leads to: 

BcosiPmD) = Ccos(pg(D - H)) 

B^sm{l3raD) = C(3gSm{(3g{D - H)) 

which implies the following dispersion relations: 

F{KlG^) = ^ta.n{/3mHT) - /3, tan(/3,i7(r - 1)) = (3) 

3. Asymptotics of the dispersion relation 

In order to give a specific example of dispersion curves, we choose the 
following parameters: h = 250nm,d = 5Qnm (i.e. r = 1/5). A colorplot of 
F{Kq,G'^) = is given in fig. 2. Two set of curves can be distinguished: 
a lower branch, below the plasma frequency of the metal (i.e. Kq < 1) and 
below the light cone {Kq — G); and an upper branch made of several curves, 
above the plasma frequency and above the light cone. The curves above the 
light cone correspond to the usual guided modes, while the mode below the 
hght cone correspond to fields localized at the air-metal interface. 

In order to be more specific, we proceed to an asymptotic study of these 
curves. Near the 'T point": (/co,7) = (0,0) and for a waveguide which is 
small with respect to Ap, that is -C 1, we get the following dispersion 
relation, from eq. (3): 

/3^/£„T-/3j(T-l)=0 (4) 

upon replacing the expression of the various quantities in terms of the prin- 
cipal variables: Kq and G, we obtain the relation: 




Figure 2: Dispersion relation of the waveguide for r = 1/5. The modes appear as yellow- 
red lines. 
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Figure 3: Dispersion relation of the waveguide for t = 1/10 and h = 50nm. The modes 
appear as yellow-red lines. The black lines correspond to the asymptotic dispersion curves 
Fu and Fd. 
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Figure 4: Dispersion relation of the waveguide for r = 1/2 and h = bOnm. The modes 
appear as yellow-red lines. The black lines correspond to the asymptotic dispersion curves 
F„ and Fd. 

This equation can be solved for Kq exactly if we assume F = 0. 
We then obtain two branches: 

Ko = + + v/4G2r + {G^ - 1)2 = (5) 

= y^^/l + G^^'^/I^^TTW^ = UG) (6) 

The limits of the branches Fu and Fd are as follows. As G tends to infin- 
ity, Fd{G) tends to the horizontal line Kq = a/1 — r. On the other hand, 
when G ~ 0, we have: Fd{G) ~ — t G. For the superior branch, when 
G ~ 0, it holds Fu{G) ~ 1 + ^G^, showing that the dispersion relation is 
parabolic, whereas when G tends to infinity: Fu{G) ~ G and the upper curve 
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Figure 5: Dispersion relation of the waveguide for r = 9/10 and h = 50nm. The modes 
appear as yellow-red lines. The black lines correspond to the asymptotic dispersion curves 
F„ and Fd. 

is asymptotic to the dispersion of the vacuum. 

In fig. (3), (4) and (5), the two branches corresponding to and are 
given, as well as the colorplot of the dispersion relation. The parameters 
are h = 50nm and r = 1/10, 1/2, 9/10. The upper branch Fu is a good 
approximation for the fundamental (non-plasmonic) mode of the waveguide, 
while the lower branch is essentially valid near G = 0. However, the filling 
ratio r = 1/2 corresponds to a situation where both the upper and lower 
curves approximate very well the true dispersion curves. 

As r varies in [0, 1], the variation of the curves is given in fig. 6. The 
part of the curve above 1 corresponds to those frequencies above the plasma 
frequency Up, in which case the real part of the permittivity of the metal 
is positive. When losses are small enough quasi-modes exist in the guide. 
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Figure 6: Bundle of the approximate dispersion curves given in (5) for a varying r between 
and 1. 

leading to the above branch of the dispersion curve. The branch is thus the 
photon-hke curve. The lower branch corresponds to the frequencies below 
the plasma frequency, for which the real part of the permittivity of the metal 
is negative. It corresponds to modes essentially localized at the interface of 
the metal with air, that is, surface plasmons. The lower branch is therefore 
the plasmon-like curve. When r tends to 0, the tangent dispersion relation 
becomes: 

Ko^l,Ko^G (7) 

and when r tends to 1: 

Ko = 0, = + 1 (8) 

In particular, the limit case r = shows that the dispersion curve at finite r 
represents an anti-crossing of the two straight lines in eq. (7), which indicates 
the strong coupling between the photon modes and the plasmon modes. 

When r = 1/2, the functions F„ and allows to obtain explicit expres- 
sions of the modes with a good approximation. In fig. 7, we have plotted 
u{x) for G = 1 and for each branches. The red curve corresponds to the 
plasmon-hke mode and the black one to the photon-hke one. It is seen that 
the red curve has the feature of a field localized at the metal-air interface. In 
fig. 8, the corresponding maps of the field ^{u{X)e^'^^) are given. The up- 
per map corresponds to the photon mode and the lower one to the plasmon 
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Figure 7: Modulus of the upper (black) and lower (red) modes at G = 1. The vertical 
blue line indicates the air-metal interface. 




Figure 8: Map of the modes for the upper and lower branches at G = 1. The modulus is 
given in fig. (7). 
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Figure 9: Dispersion curves for r = 0.05 and h = 200nm. The curvature of the plasmonic 
branch can be seen. 

mode. The concentration of the field near the metal-air interface is clearly 
seen. Finally, an interesting situation arises for small filling ratios. For in- 
stance for r = 1/10, it can be seen in fig. (3) that the plasmonic branch 
is curved, in such a way that for some value of Ko (near 0.8) there are two 
corresponding values of G satisfying the dispersion relation. Therefore by 
exciting both modes at their eigenfrequency, one obtains a spatial oscillation 
between these two modes. Indeed, denoting ipi and ip2 these modes, at the 
corresponding frequency, a general mode is a linear combination of ipi and 
-02 in the form: 

which in modulus takes the usual form of the beating phenomenon: 

Z)\^ = |aV'iWr + |/3V'2(X)|2 + 2a/3^i(X)^2(X) cos ((d - ^2)^) 

Specifically, we choose Kq = 0.72, solve for F{Kq,G'^) = and find Gi ~ 
0.9556 , G2 ~ 2.8563. The modulus of the modes is given in fig. (10) and the 
map of the modes are given in fig. (11). It can be seen on these figures that 
the smaller spatial period is also much more concentrated at the air-metal 
interface than the other one. This shows that, despite the fact that the layer 
of metal is very thin and that the metal is near the plasma frequency, it is still 
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Figure 10: Modulus of the two plasmonic modes obtained a.t K = 0.72 of. fig. (9). 
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Figure 11: Map of the two plasmonic modes whose modulus is given in fig. (10). 
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possible to obtain plasmon-like modes strongly localized near the air-metal 
interface. The fact that for a small enough filling ratio r the plasmonic 
branch is not monotonic with G, i.e. has a maximum and bends is by itself 
an interesting phenomenon that seems to indicate that the dispersion band 
can be described by an effective, spatially dispersive, medium. It should be 
noted that the extreme thinness of the layer of metal questions the validity 
of the very simple model used here. In fact, such a thin layer is strongly 
corrugated and the use of a Drudc model is questionnable. A more complex 
model is needed, for instance by describing the surface of the metal as a 
random rough surface [8, 9, 10], to draw sound conclusions. Work is in 
progress in that direction. 

4. Some remarks on the modes 

When there are no losses in the layer of metal, the spectral problem 
defined by equation (1) provides a complete set of eigenvectors parametrized 
by G. When r = (i.e. no metal), this corresponds to the following set: 

i^niX) = ^J^cos{^),Ko = ^^2 + ^ with the inner product: {f,g) = 

J^'' f{X)g{X)dX. 

The adjunction of losses implies that the solutions to F{Kq, G"^) = are 
no longer real. However, it is still possible to choose a real Kq and obtain 
a complex G, or, reciproquely, to choose a real G and obtain a complex Kq 
(cf. the discussion in [11]). To that extend, the true situation is that one can 
choose a complex Kq or a complex G and obtain a corresponding G or Kq. 
In fact, the relation F{Kq, G"^) = defines a multi-sheeted Riemann surface. 
An explicit example is given in fig. 12. 

The modes read as: 

H(X, Z,t) = 3fJ(H(X)e'(^^-^°^^)) 

Therefore, whenever G or Kq is complex, this leads to an exponential damp- 
ing in one direction of time or space and an exponential increasing in the 
other direction. These modes are no longer tempered distribution in both 
time and space. In the usual hilbertian framework of square integrable func- 
tions L^{[Q,H] X M), there are no longer modes extending to infinity in an 
infinite lossy structure. A more convenient mathematical setting would be 
that of hilbertian triad, or Gel'fand triple [12, chap. 1]. Still, the matching of 
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Figure 12: Four-sheeted Riemann surface corresponding to the approximate dispersion 
relations given in eq. (5). 

two such modes can be bounded. Consider as a very simple example the func- 
tions: and e^. By matching them at x = 0, we obtain the function e~'^'. 
This what happens with Laplace transform, which provides a decomposition 
over a half space: f{x) = J^°° F{p)e~^^dx which is convergent whenever 
3fJ(p) > 0. We can apply the same kind representation in the lossy waveg- 
uide. Indeed, in such a structure the sources of the fiel cannot be placed 
at infinity as explained above. The sources here are of two kind: localized 
sources inside the waveguide such as quantum dots, atoms... and distributed 
sources inside the metal accounting for the fluctuations due to the absorp- 
tion. These last can be accounted for by introducing Langevin current J such 
as in [13]. The other can be obtained by using sums of Green function: the 
localized sources can be modelized as Dirac delta functions. When only the 
X dependance is considered an infinite set of energies are obtained at 
fixed G. When losses are introduced, the modes can be written in the form 
ipni^', r) where ipni^', 0) are the modes of the lossless guide. For a fixed X, 
these functions are all analytic in the variable F, which implies that, as the 
family is complete at F = 0, it remains complete for a small enough F [14]. In 
a forthcoming work, we shall give explicit expressions for the Green function 
and apply it to the cavity quantum electrodynamics of plasmons. 
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5. Conclusion 



We have studied the effect of confinement on the surface plasmons that 
can exist at an air-metal interface. We have given some approximate explicit 
solutions to the dispersion curves which has allowed us to demonstrate the 
strong coupling between the plasmon-like and the photon-like modes. We 
have shown the existence of strongly localized plasmon-like modes and also 
the onset of spatial dispersion when the filling ratio of metal is very low. 
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